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Abstract 

The ^He beta decay and ''Be electron capture processes are studied using 
variational Monte Carlo wave functions, derived from a realistic Hamiltonian 
consisting of the Argonne fig two-nucleon and Urbana-IX three-nucleon inter- 
actions. The model for the nuclear weak axial current includes one- and two- 
body operators with the strength of the leading two-body term — associated 
with A-isobar excitation of the nucleon — adjusted to reproduce the Gamow- 
Teller matrix element in tritium /3-decay. The measured half- life of ^He is 
under-predicted by theory by ~ 8%, while that of ^Be for decay into the 
ground and first excited states of ^Li is over-predicted by ~ 9%. However, 
the experimentally known branching ratio for these latter processes is in good 
agreement with the calculated value. Two-body axial current contributions 
lead to a ~ 1.7% (4.4%) increase in the value of the Gamow- Teller matrix ele- 
ment of ^He (''Be), obtained with one-body currents only, and slightly worsen 
(appreciably improve) the agreement between the calculated and measured 
half- life. Corrections due to retardation effects associated with the finite lep- 
ton momentum transfers involved in the decays, as well as contributions of 
suppressed transitions induced by the weak vector charge and axial current 
operators, have also been calculated and found to be negligible. 

23.40.Bw, 23.40.Hc 
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I. INTRODUCTION 



The present study deals with weak transitions in the A=Q and 7 systems within the 
context of a fully microscopic approach to nuclear structure and dynamics, in which nucleons 
interact among themselves via realistic two- and three-body potentials, and with external 
electroweak probes via realistic currents consisting of one- and many-body components. To 
the best of our knowledge, calculations of the superallowed ^He beta-(/3-) decay and ^Be 
electron- (e-) capture processes have relied in the past on relatively simple shell- model or 
two- and three-cluster wave functions. The shell model calculations have typically failed to 
reproduce the measured Gamow- Teller matrix elements governing these weak transitions, 
unless use was made of an effective one-body Gamow- Teller operator, in which the nucleon's 
axial coupling constant qa is quenched with respect to its free value — for a recent summary 
of a shell-model analysis of /5-decay rates in A < 18 nuclei, see Ref. 

More phenomenologically successful models have been based on aNN (for A=Q) or a-t 
and a-T (for A=7) clusterization, and have used either Faddeev techniques with a separa- 
ble representation of the NN and aN potentials or the resonating-group method [Q. 
However, while these studies do provide useful insights into the structure of the A=6 and 
7 nuclei, their connection with the underlying two- (and three-) nucleon dynamics is rather 
tenuous. For example, it is unclear whether the required quenching of qa in the shell-model 
calculations reflects deficiencies in the associated wave functions — a lack of correlations or 
limitations in the model space — and/or in the model for the axial current operator — which 
typically ignores many-body terms — or, rather, a true modification of the nucleon axial 
coupling in medium. 

In this work, we use variational Monte Carlo (VMC) wave functions 0-0] ? derived from 
a realistic Hamiltonian consisting of the Argonne fig two-nucleon and Urbana-IX three- 
nucleon ^ interactions. The VMC wave functions provide a reasonable description of the 
energy spectra of low-lying states in A=6-8 nuclei P,p!0| , and of elastic and inelastic elec- 



tromagnetic form factors and radiative widths of Li states ||TT 



The model for the nuclear weak vector and axial- vector currents is that developed in Refs. 
r^,|TB|, consisting of one- and two-body terms. The weak vector charge and current operators 



are constructed from their isovector electromagnetic counter-parts |TJ,|T5|, in accordance 
with the conserved-vector-current (CVC) hypothesis. The leading two-body term in the 
axial current is due to A-isobar excitation, while the leading two-body axial charge operator 



is associated with the long-range pion-exchange term [16|, required by low-energy theorems 



and the partially-conserved-axial-current relation. The largest model dependence is in the 
weak axial current. To minimize it, the N ^ A transition axial coupling constant has been 
adjusted to reproduce the experimental value of the Gamow- Teller matrix element in tritium 
/?-decay in an essentially exact calculation, using correlated-hyperspherical-harmonics wave 
functions, derived from the same Hamiltonian adopted here . 

This manuscript falls into seven sections. In Sec. |I| explicit expressions for the /5-decay 
and e-capture rates are derived in terms of reduced matrix elements of multipole operators 
of the nuclear weak current, whose model is succinctly described in Sec. 0. The A=6 
and 7 nuclei VMC wave functions are reviewed in Sec. |T|, while predictions for the ®He 
/3-decay and ''Be e-capture rates are presented and discussed in Sec. [V^. Our conclusions 
are summarized in Sec. [VII . 
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II. TRANSITION RATES 



Nuclear /3-decay is induced by the weak-interaction Hamiltonian 

G 
72 



i7;^ = ^|dxe-(P^+P^)-^/.nx), 



(2.1) 



where Gy is the Fermi couphng constant (Gy =1.14939 10 ^ GeV ^ 10), ^o- is the leptonic 
weak current 



(2.2) 



and j"{x.) is the hadronic weak current density. The electron and (electron) antineutrino 
momenta and spinors are denoted, respectively, by pe and p^, and Ug and (e~ emission 
is being discussed here). The Bjorken and Drell conventions are used for the metric 
tensor g'^'^ and 7-matrices. However, the spinors are normalized as ulu^ = v^Vy = 1. The 
amplitude for the process "^Z {Z + l)e~Ve is then given by 



{mw\^) 



G 



V 



V2 



r(-q;^(Z + 1), Jf Mf\jl{q)\^Z, J, M,) 



(2.3) 



where (i = pe + Pu, \^Z, Ji Mi) and \^{Z + 1), ^/ Mf) represent the initial and final nuclear 
states, the latter recoiling with momentum — q, with spins Ji and Jj and spin projections 
Mi and Mf, respectively, and 



r(q) = / dxe^''-j-(x) = (p(q),j(q)) . 



(2.4) 



Standard techniques [|15,|1^ are now used to carry out the multipole expansion of the 
weak charge, p(q), and current, j(q), operators in a reference frame in which the z-axis 
defines the spin-quantization axis, and the direction q is specified by the angles 6 and (p: 

(2.5) 
(2.6) 

{JfMf\e;^,.i\cO\J.M,) = -V2^Y.Ki [±Miiq) + Eiiq)] , (2.7) 

i>i 



{JfMf\p\q)\JiMi) = v^^X^a(g) 



l>0 



Jf Mf I e^o ■ jt (q) I J, M,) = ^ (g) 



l>0 



where the dependence on the direction q, and on the initial and final spins and spin projec- 



tions is contained in the functions X[, with A=0, ±1, defined as 



21 + 1^ 



1/2 



Dl,i-(l),-e,0){JiM,,lQJfMf) , (2.8) 



with lz=Mf — Mi. In Eqs. ( p.5|) - (|2.7|) the g-dependent coefficients are the reduced matrix 
elements (RME's) of the Coulomb (C), longitudinal (L), transverse electric {E) and trans- 
verse magnetic (M) multipole operators, as given in Refs. ||13|J18[| , and the vectors b\ denote 
the linear combinations 
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eg3 , 
1 



(2.9) 
(2.10) 



with 6^3 = q, 6^2 = z X q/|z x q|, e^i = x 6^3. Since the weak charge and current 
operators have scalar/polar- vector {V) and pseudoscalar/axial- vector {A) components, each 
muhipole consists of the sum of V and A terms: 



Ta^=TuAV)+TaM) 



(2.11) 



where T=C, L, E, and M, and the parity of /th-pole \^-operators is opposite of that of Zth- 
pole A-operators. The parity of Coulomb, longitudinal, and electric /th-pole V^-operators is 
(— )', while that of magnetic /th-pole V^-operators is (— )'+^. Finally, in Eq. ( p.8|) the Dl^ 
are rotation matrices in the standard notation of Ref . |^ . 
The rate for nuclear beta decay is then obtained from 



(2,12) 



where Ei=mi is the rest mass of the initial nucleus '^Z, Ef=py + \JpI + + ^Jcf + mj 
is the energy of the final state, rUe and ruf being the rest masses of the electron and final 
nucleus "^{Z + 1). Carrying out the spin sums leads to ||13|,|18| 



^2X + 1 



;i + v,-v.)^|Q(g)P 

l>0 

+ (1 - Ve-v^ + 2Ve-qv^-q)^|Li(g)|2 

l>0 

-2q.(ve + v,) Y.^e[Ci{q)L*{q)] 

l>0 

+ (1 - v,-qv,-q) 5: [|MKg)P + 
i>i 

-2q-(ve-v,)5:Re[M,(g)E;(g)] 



i>i 



(2.13) 



where Ve=Pe/yPe + ml and Vj,=p,^ are the velocities of the electron and neutrino, respec- 
tively. Note that angular momentum and parity selection rules restrict the number of non- 
vanishing /-multipoles contributing to the transition. 

A few comments are now in order. Firstly, the expression above is valid for e~ emission, 
however, the rate for emission corresponding to a transition ^Z — > ^{Z — l)e^z/e has 
precisely the same form, but for the sign in the last term of Eq. ( p.l3| ), +2 q-(ve — Vj,) ■ ■ ■ . 
Secondly, for allowed transitions with \ Ji — Jf \ = ±1, and vTj TTf=l the above rate is easily 
shown to reduce to the familiar expression in terms of Fermi [F (x Co{q = 0; V)] and Gamow- 
Teller [GT oc Ei(q = 0; A)=^/2 Li{q = 0] A)] matrix elements, see Sec. [V| and Ref 



m 



Finally, a more realistic treatment — like that outlined in Ref. |]2^, for instance — takes into 
account the distortion of the outgoing e"^ wave function in the Coulomb field of the residual 



4 



nuclear system. The simple approximation of multiplying the right-hand-side of Eq. (p.l3|) 



by the ratio of the charged lepton density at the nuclear radius to the density at infinity has 
been deemed, however, to suffice for our purposes here (see Sec. [V]). One of the objectives 
of the present study is to estimate, in an allowed decay such as the ^He — Lie~I7e under 
consideration [with ( Jf^ Ti) = (0+;1) and ( Jj^; T/) = (l+;0)], the size of corrections associated 
with i) retardation effects due to the finite lepton momentum transfer in the decay, and ii) 
transitions other than those of F [Co(g; V^)] and/or GT [Li{q; A) and Ei{q] A)] type. For 
example, a naive analysis of the 6-body decay above would ignore the momentum transfer 
dependence — it is in the range < g < 4 MeV/c — as well as the contributions arising 
from transitions induced by the axial charge and vector current operators via Ci(g; A) and 
Mi{q; V) RME's, respectively. Of course, available tabulations of F and GT matrix elements 
extracted from experiment do attempt to estimate these corrections (as well as those 
due to electron screening, the finite extent of the nuclear charge distribution, etc.). The 
latter are typically obtained, however, within the context of a shell-model description of 
the nuclear wave functions. Thus, it is interesting to re-examine the issue above within the 
present framework using more realistic wave functions. 

The e~ capture (so-called e-capture) process is governed by the same weak-interaction 
Hamiltonian in Eq. ( [2.11 ). However, the lepton weak-current density is now given by 

Z,(x) = e-^P''-''n^7,(l - 75) (Atom/ | ^e(x) | Atom^) , (2.14) 

where ipe{^) is the electron field operator, and | Atom^) and | Atom/) are the initial and 
final atomic states, respectively. A realistic description of nuclear electron capture requires. 



therefore, a careful treatment of the atomic physics aspects of the process ||2^, such as those 
relating to atomic wave-function overlaps, exchange contributions, and electron-correlation 
effects. We defer to Sec. |V| for a discussion of some of these issues in the context of the ^Be 
e-capture of interest here. In this section, however, we simply approximate 

(Atom/ I ^e(x) I Atom,) ^ ^^^x(^e) = ^^u{p,, s,) ^ , (2.15) 



ignoring atomic many-body effects and relativistic corrections — this latter approximation is 
reasonably justified for low Z, since the e~ velocity is ~ Za ^ 1. In Eq. ( p.l5| ), Ris{x) is 



the Is radial solution of the Schrodinger equation, and the two-component spin state xi^e) 
has been conveniently replaced by the four- component spinor u{pe, Sg) in the limit Pe — 0, 
which allows us to use standard techniques to carry out the spin sum over Sg at a later stage. 
The transition amplitude reads 

{f\Hw\z) = -^i?i,(0)r(-p,;^(Z-l),J/M/|jt(p,)|^Z,J,M,) , (2.16) 



5vr 

where la = Uu'jai^ ~ l5)ue, and the e~ radial wave function has been factored out from the 
matrix element of j'^{pu) by approximating it with its value at the origin. The resulting 
differential rate is then written as [0,0] 

dF}, = 2vr5(E, - Ej) ^ ^ | {f\H^\t) \' |^ , (2.17) 



MiMf ■ 
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where the initial and final energies are now given by Ei=mi + nie and Ef=py + yp^ _|_ 
respectively. Note that atomic binding energy contributions have been neglected, since they 
are of the order (Za)^me/2. The square of the amplitude summed over the spins can be 
obtained, mutatis mutandis, from that corresponding to e"*" emission discussed above in the 
limit pe=0 (q=p^): 



Ji -^ L) g^g^ 



An 



An (2 J, + 1) 



l>0 



+ Y.\Mi{Pu)-Ei{p,) 



i>i 



(2.18) 



III. WAVE FUNCTIONS 



The VMC wave function, "^t^J^', T), for a given nucleus, is constructed from products of 
two- and three-body correlation operators acting on an antisymmetric single-particle state 
with the appropriate quantum numbers. The correlation operators are designed to refiect 
the infiuence of the interactions at short distances, while appropriate boundary conditions 
are imposed at long range f^-^]- The \E't(<^'^; 7") has embedded variational parameters that 
are adjusted to minimize the expectation value 



(vE^tI^t) 



> En 



(3.1) 



which is evaluated by Metropolis Monte Carlo integration. 
A good variational trial function has the form 



1^1 



1+ E 

i<j<k 



i<j 



(3.2) 



The Jastrow wave function, \E'j, is fully antisymmetric and has the {J'^]T) quantum numbers 
of the state of interest. For the s-shell nuclei we use the simple form 



II fijk 
i<j<k 



II Mr 



\<^a{JMTTs)) 



(3.3) 



Here fd^ij) and /^^-^ are central two- and three-body correlation functions and $^ is a Slater 
determinant in spin-isospin space, e.g., for the a-particle. 



|$4(0000)) = Alp'i p In^ n I) . 



(3.4) 



The Uij and U/j^ are noncommuting two- and three-nucleon correlation operators, and S 
indicates a symmetric sum over all possible orders. The Uij includes spin, isospin, and tensor 
terms: 



Uij = E Mr^j)Ofj 
p=2,6 



(3.5) 
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where the Of~^'^ = [1, Uj ■ aj, Sij] [1, • tj] are the main static operators that appear in the 
NN potential. The fc{r) and Up{r) functions are generated by the solution of a set of coupled 
differential equations which contain a number of variational parameters The has 
the spin-isospin structure of the dominant parts of the NNN interaction as suggested by 
perturbation theory. 

The optimal Uij and U™^ do not change significantly from nucleus to nucleus, but 
"^j does. For the p-shell nuclei, \Ef j includes a one-body part that consists of 4 nucleons 
in an a-like core and [A — 4) nucleons in p-shell orbitals. We use LS coupling to obtain 
the desired JM value, as suggested by standard shell- model studies [|^. We also need to 
sum over different spatial symmetries [n] of the angular momentum coupling of the p-shell 
nucleons The one-body parts are multiplied by products of central pair and triplet 

correlation functions, which depend upon the shells {s or p) occupied by the particles and 
on the LS[n] coupling: 

i^i)=^{ n f^k n fss{r^,) n /^^m 

i<j<k «<j<4 k<4:<l<A 

E {f^LSM n fif^''\rirn)\^AiLS[n]JMTTs)l23A-.,...A)) } . (3.6) 

LS[n] 4:<l<m<A 

The operator A indicates an antisymmetric sum over all possible partitions into 4 s-shell 
and {A — 4) p-shell particles. The pair correlation for both particles within the s-shell, fss, 
is set to the fc of the a-particle. The fgp is similar to fss at short range, but with a long- 
range tail going to a constant ~ 1, while the /j^'^'"^ is allowed to depend on the particular 
single-particle channel. 

The LS[n] components of the single-particle wave function are given by: 

|<I>a(L5MJMTT3)i234:5...a) = |<f4(0000)i234) H (3-7) 

4<1<A 

{ [ n Ylmi^al)]LM^[n] X [ [] Xli^ms)]sMs \ X [ [] Z^K|^3)]Tr3) • 
[ 4<1<A 4<1<A ) jj^j 4<1<A 

The (j)p^{Rai) are p-wave solutions of a particle in an effective a-N potential that has Woods- 
Saxon and Coulomb parts. They are functions of the distance between the center of mass 
of the a core and nucleon /, and may vary with LS[n]. The wave function is translationally 
invariant, so there is no spurious center of mass motion. 

Two different types of \E'j have been constructed in recent VMC calculations of light 
p-shell nuclei: an original shell-model kind of trial function P which we will call Type I, 
and a cluster-cluster kind of trial function which we will call Type II. In Type I trial 
functions, the (pp^ir) has an exponential decay at long range, with the depth, range, and 
surface thickness of the Woods-Saxon potential serving as variational parameters. The fsp 
goes to a constant somewhat less than unity, while fpp^"^^ is similar to fss at short-range 
with an added long-range tail that is larger for states of lesser spatial symmetry [n]. Details 
for these A = Q,7 trial functions are given in Ref. [Q. 

In Type II trial functions, 4>p^{r) is again the solution of a p-wave differential equation 
with a potential containing Woods-Saxon and Coulomb terms, but with an added Lagrange 
multiplier that turns on at long range. This Lagrange multiplier imposes the boundary 
condition 
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^(r ^ oo)]" (X l^fc„(27r)/r, (3.8) 

where WkmC^lf) is the Whittaker function for bound-state wave functions in a Coulomb 
potential and n is the number of p-shell nucleons. The 7 is related to the cluster separation 
energy which is taken from experiment. The accompanying fgp goes exactly to unity (more 
rapidly than in the Type I trial function) and the fpp^^^ are taken from the exact deuteron 
wave function in the case of ^Li, or the VMC triton (^He) trial function in the case of ^Li 
(''Be). Consequently, the Type II trial function factorizes at large cluster separations as 

^aV'rW^fcm(27rar)/rar . (3.9) 

where ipa and ipr are the wave functions of the clusters and r^^ is the separation between 
them. More details on these wave functions are given in Refs. 0,01 • In the case of ^He, 
which does not have an asymptotic two-cluster threshold, we generate a fpp^""^ correlation 
assuming a weakly bound ^Sq nn pair. 

For either type of trial function, a diagonalization is carried out in the one-body basis 
to find the optimal values of the Pisin] mixing parameters for a given {J^]T) state. The 
trial function, Eq. (|3.2| ), has the advantage of being efficient to evaluate while including the 
bulk of the correlation effects. A more sophisticated variational function can be constructed 
by including two-body spin-orbit correlations and additional three-body correlations P,p7 



but the time to compute these extra terms is significant, while the gain in the variational 
energy is relatively small. In calculations of these extra terms have a negligible (< 0.2%) 
effect on the GT matrix element. 

The wave function at a given spatial configuration R = ri,r2, ...ta can be represented 
by a vector of 2^ x I {A, T) complex coefficients in spin and isospin space [^]. For the nuclei 
considered here, this gives vectors of 320 in ^Li, 576 in ^He, and 1,792 in ^Li and ^Be. The 
spin, isospin, and tensor operators 0^~^'^ contained in the Hamiltonian and other operators 
of interest are sparse matrices in this basis. 



IV. NUCLEAR WEAK CURRENT 



The model for the nuclear weak current has been most recently and exhaustively de- 
scribed in Ref. |T^. Here we summarize only its main features. 

The nuclear weak current consists of vector and axial-vector parts. 



p±(q) =p±(q; ^) +p±(q;^) 
j±(q) =j±(q;^) +j±(q;^) , 



(4.1) 
(4.2) 



with corresponding one- and two-body components. The weak vector current is constructed 
from the isovector part of the electromagnetic current, in accordance with the conserved- 
vector-current (CVC) hypothesis 



j±(q;'^)= r±Jz(q;7) 



(4.3) 



where j±(q; V) is the charge-lowering (-) or charge-raising (+) weak vector current, izi^l'i l) 
is the isovector part of the electromagnetic current, and T± is the (total) isospin- lowering 
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or isospin-raising operator. A similar relation holds between the electromagnetic charge 
operator and its weak vector counterpart. For reference, we list only the expressions for the 
one-body terms in p(q; V) and j(q; V), in the notation of Ref. [|1^]: 

pS'Hq; V) = pi^M. V) + p!Sc(q; v) , (4.4) 

with 

pSrIq; V) = e^'^-^^ , (4.5) 



pSc(q;^) = -i 



and 



.(2/i^-l) 



4^2 



Pi , e 



i"i,± q • (cTi X Pi) e 



— i—Ti ± q X cTj e'''''"* 



(4.6) 



(4.7) 



where [■ ■ ■ , ■ ■ ■]+ denotes the anticommutator, p, cr, and r are the nucleon's momentum, 
Pauh spin and isospin operators, respectively, and fi^ is the isovector nucleon magnetic 
moment (/i^ = 4.709 hn)- Finally, the isospin raising and lowering operators are defined as 



(4.8) 



The one-body terms in the axial charge and current operators have the standard ex- 
pressions ||13[ obtained from the non-relativistic reduction of the covariant single-nucleon 
current, and retain terms proportional to m being the nucleon mass: 



and 



with 



A) = i\:>M, A) + i\:ic{ci; A) 



9a 
2 m 



(1) 



j!,NR(q; ^) = -9a ri,± (Ti e'^i-' 



(4.9) 



(4.10) 



(4.11) 



•(1) , 



q;A) 



9A 

1 

-2^ 



n,± ( 


P- , e^'^-^' 




(Ji ■ Pi Pi , e"' 






+ 





1 



■ q 



Pi , e 



iq-fi 



o-i ■ Pi , e 



iqri 



+ i q X Pi e 



iqi-i 



iqr, 

Ti,±qcri ■ qe'^ * 



2mm 



(4.12) 



The axial-vector coupling constant qa is taken to be |^| 1.2654±0.0042, by averaging 
values obtained from the beta asymmetry in the decay of polarized neutrons and the half- 
lives of the neutron and superallowed 0"^ 0"*" transitions. The last term in Eq. ( |4.12| ) is the 
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(m^ is the muon mass) 



for which the couphng constant 



induced pseudoscalar contribution 
gp is taken as gp=-6.78 qa- 

Some of the two-body axial-current operators are derived from vr- and p-meson exchanges 
and the pvr-transition mechanism. These mesonic operators, first obtained in a system- 
atic way in Ref. |3y], have been found to give rather small contributions in weak transi- 



tions involving few-nucleon systems [lT2| , p^ , p!7| . The two-body weak axial-charge operator 
includes a pion-range term, which follows from soft-pion theorem and current algebra argu- 
ments [p!6| , pT| , and short-range terms, associated with scalar- and vector-meson exchanges. 
The latter are obtained consistently with the two-nucleon interaction model, following a 
procedure similar to that used to derive the corresponding weak vector-current opera- 
tors n. 

The dominant two-body axial current operator, however, is that due to A-isobar ex- 
citation ||T2| , [l3[| . Since the A^A transition axial- vector coupling constant g\ is not known 
experimentally, the associated contribution suffers from a large model dependence. To re- 
duce it |ll2| , pr3| , the coupling constant g\ has been adjusted to reproduce the experimental 
value of the Gamow- Teller (GT) matrix element in tritium (3 decay, 0.957 ± 0.003 |0. The 
value used for g\ in the present work is g\- 



--\.VJ gA 

The A-excitation operator used here is that derived, in the static A approximation, using 
first-order perturbation theory. This approach is considerably simpler than that adopted 
in Ref. |jl3|, where the A degrees of freedom were treated non-perturbatively, by retaining 
them explicitly in the nuclear wave functions [Q. However, it is important to emphasize 
that results obtained within the perturbative and non-perturbative schemes are within 1% of 
each other typically, once g\ is fixed, independently in the perturbative and non-perturbative 
calculations, to reproduce the experimentally known GT matrix element; see Table XV in 
Ref. m 



V. CALCULATION 

The calculation of the /3-decay and e-capture rates proceeds in two steps: firstly, the 
Monte Carlo evaluation of the weak charge and current operator matrix elements, and the 
subsequent decomposition of these in terms of reduced matrix elements (RME's); secondly, 
the evaluation of the rate by carrying out the integrations in Eqs. (|2.12|) and ( p.l7|) . 

The RME's are obtained from Eqs. ( |2.5| )- P^ ) by choosing appropriately the q-direction. 
For example, in the /3-decay of ^He the Ci(g; A) and L\[(i] A) RME's are determined from 



CM: A) 



fLi,10|pU9z;A)|6He,00) , 
fLi,10|z-jUgz; A)|6He,00) 



while the Ei {q; A) and Mi (g; ) RME's are determined from 

E^{q■ A) = -^('Li, 10|z ■ jt(gx; A)f}le, 00) 



M^{q■ V) 



1 



{'U,10\y-il{q±-V)\^}ie,00) . 



(5.1) 
(5.2) 

(5.3) 
(5.4) 
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Here Ji,Mj=0,0 and J/,Mj=l,0 and the spin-quantization axis is along z. The matrix 
elements above are computed, without any approximation, by Monte Carlo integration. 
The wave functions are written as vectors in the spin-isospin space of the A nucleons (^4=6 
or 7 here) for any given spatial configuration R = (ri, . . . , ta). For the given R, the state 
vector 0^(R)\['j(R) is calculated with techniques similar to those developed in Ref. [|1^ — 
0(R) is any of the operators p(q; V), j(q; V), etc., and \l/i(R) is the wave function of the 
initial nucleus. The spatial integrations are carried out with the Monte Carlo method by 
sampling R configurations according to the Metropolis algorithm, using a probability density 
proportional to (\l'|(R)\E'j(R)), where \1//(R) is the wave function of the final nucleus and 
the notation (■ ■ ■) implies sums over the spin-isospin states. Typically 20,000 configurations 
are enough to achieve a relative error < 1 % on the matrix elements. 

Once the RME's have been obtained, the calculation of the total transition rate is reduced 
to performing the integrations over the outgoing momenta in Eqs. (|2.12| ) and (|2.17| ). A glance 
at Eqs. ( |2.12| ) and ( ^.l^j| ) shows that the differential rate depends on the magnitude of the 
electron momentum and the variable Xeu='Pe ■ Pu, since the magnitude of the neutrino 
momentum is fixed by the energy- conserving 5-function. The total rate can therefore be 
written as 



Gy 



fi 



2 7r3 



(5.5) 



where the constant is simply given by 



/3 



27r 



2Ji + l Jo 



(5.6) 



Here all momenta and rest masses have been expressed in units of the electron mass 



{Pe=Pe/me 



rrii 



--rrii/me 



TUf, 



-1, etc.^ 



denotes the content of the large square brackets 
on the right-hand-side of Eq. ( p.l3|) , is the recoil factor resulting from integrating out 
the ^-function, 



/rec " 

and the neutrino momentum p^^ is 



^ Pe -^eu Pu 

rfif rfif 



2 A 



(5.7) 



1 +p^Xeu/mf + y (1 + p^Xeu/mfY + 2 A/m/ 



(5.8) 



where A=mi — rfif — \jp\ + 1 —pl/{2mf). Lastly the function F{Zf,p^) accounts approxi- 
mately for wave-function distortion effects of the outgoing electron in the Coulomb field of 
the final nucleus with atomic number Zf and radius Rf (expressed in units of the Compton 
wavelength of the electron) |^ 



F(Z,pJ = 2(l + 7o) i2p,R) 



Tn m2(70-l) 



|r(7o + iu) 
|r(27o+l) 



(5.9) 
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with 7o = Y 1 — {aZy and u = aZ/ve- The maximum allowed electron momentum is 
denoted with p*, the upper integration limit in Eq. (|5.6| ), and is given by 



— * 

Pe 



mj + 2mf {rrii — rrif) + 1 — ruf — 1 . (5.10) 



2 



Neglecting the recoil of the final nucleus leads to: /rec=l, Py='rrii — nif — ^pl + 1, and 
Te = \J{'mi-mff - 1. 

Note that there is an implicit dependence on pe and Xev in the RME's via the momentum 
transfer q. It is convenient to make this dependence explicit by expanding the RME's as 

Tz(g)=g'"E^/,2ng'", (5.11) 

n>0 

consistently with the known expansions of the multipole operators in powers of q |]13[ . Here 
m = Z, Z ± 1, depending on the RME considered. For example, in the ^He /3-decay one has 
Ci(g; A)=q{cifi + Ci^2(f + ■ ■ ■)' Li[q; A)=lifi + /i,25'^ + ■ ■ etc. . Given the low momentum 
transfers involved in all transitions under consideration, g < 4 MeV/c, the leading and next- 
to- leading order terms ti^ and are sufficient to reproduce accurately T/(g). Incidentally, 
the long-wavelength-approximation corresponds, typically, to retaining only the tifl term. 
Finally, standard numerical techniques — Gaussian quadratures — are used to carry out the 
integrations in Eq. ( ^.6|) . 

It is useful to consider the case of allowed transitions for which |Jj — J/|=±1,0 and 
7rj7rj=l, such as the ^He /3-decay of interest in the present study. Ignoring retardation 
corrections due to the finite momentum transfer involved in the decay, one finds that the 
only surviving RME's in the limit g — *• are CoiV), Li{A), and Ei{A) — of course, CoiV) 
vanishes unless Ji=Jf. The one-body terms in the associated multipole operators read in 



this limit 13,18 



Coo{q = 0;V) = ^ (5-12) 



E,i^{q = 0;A) = V2L,i^{q = 0-A) = -^gAY^^,±^iAh , (5-13) 

VOTT 



where ai^u^ denote the spherical components of cTj. The constant 7^. is then obtained as, 
neglecting nuclear recoil corrections, 

, \Fl + g\jGTl 

2J. + 1 ^ ^'-'^^ 



rm^-mf , „ \ 1/2 / \2 

f = J^ dE,E,[E^-l) {m,-mf-E,) F{Zf,E,) , (5.15) 

since the integration over x^u can now be performed trivially. The familiar definitions of the 
Fermi and Gamow- Teller (reduced) matrix elements. 
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i 

GT = {Jf\\Y,%±o-i\\Ji) 



(5.16) 
(5.17) 



have been introduced, as well as Ee=ypl + 1. Combining Eqs. (|0| ) and (|5.14|) leads to the 
standard expression of the decay rate for allowed transitions as, for example, in Ref. p3 |. 
Finally, the total rate for e-capture easily follows from Eqs. (|2.171 ) and (|2.18|) : 



^2 ^5 



—2 



fi 



4-^-2 



E:\Ru{0)rfn 



An 



2Ji + l 



(5.18) 



where now ■ ■ ■ denotes the content of the large square brackets on the right-hand- 
side of Eq. ( |2.18| ), and the recoil factor is given by f~^=\l — Ey/iWii + 1)|~ 1 with 



E 



{rrii + 1) — m 



I [2(mj + 1)] ^ mj + 1 — rrif. Again, we have expressed masses and 

energies in units of me, and i?is(0)=mg/^i?is(0). In particular, for allowed transitions and 
ignoring recoil effects, we obtain the familiar result 



Gv K I o 



4 7r2 



EARism 



\F\^ + g\\GT\ 
2 J, + 1 



(5.19) 



As already pointed out in Sec. 0, however, a more accurate treatment of the atomic 
physics aspects is warranted for a meaningful comparison with experiment. In the case of 
the ^Be e-capture of interest here, such a program has indeed been carried out by Chen and 
Crasemann |34|. They use multi-configurational Hartree-Fock (MCHF) wave functions to 



represent the initial ground state of the 4Be atom as 

= Ci $(ls2 2s^) + C2 $(1^2 2p2) 



(5.20) 



and the final Is- or 2s-hole states, after K or L/ captures in standard nomenclature p4| 



as 



= C[ $'(15 2s2) + C'^ $'(is 2/) 
= $'(ls^2s) . 



(5.21) 
(5.22) 



They then proceed to evaluate the matrix elements (\E'i^|'?/'e(x=0)|\E'j) and similarly for 
"^L. The MCHF approach goes beyond the independent-particle approximation, commonly 
adopted in the analysis of electron capture, by retaining the effects of Coulomb correlations 
among the electrons. The end-result is that the rate, including both K and Lj captures, can 
be conveniently written as in Eq. ( ^.18| ), but for the replacement |i?is(0)p ^ B x \Ris 
where 



B = Bk + 



i?i.(0) 



Bt 



(5.23) 



in the notation of Ref. |34|, with Ba=\{^ a\'ipe{'^=0)\'^i)\'^ and a=K or L. Using the values 
listed in Table III of Ref. for B^, B^, and the ratio of (radial) wave functions at 
the origin, as well as the value for -Ris(O) from Table IX of Ref. ||2^, we find the relevant 
combination B x \Ris{0)\'^ to be equal to 7.2403 x 10~^ for the case of '''Be. 
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VI. RESULTS 



In this section we report results for the /3-decay of ^He and the e-capture in ''Be. The 
variational Monte Carlo (VMC) wave functions of the A=6 and 7 nuclei have been obtained 
from a realistic Hamiltonian consisting of the Argonne fig two-nucleon [§] and Urbana-IX 
three-nucleon interactions, the AV18/UIX model. 

The binding energies and radii predicted by the VMC wave functions of Type I and II 
discussed in Sec. |T| are listed for reference in Table |, along with the exact GFMC energies 
for this Hamiltonian and the experimental values. We note that the Type I wave function 
gives a slightly better energy in all cases, except for ^Li, while the nucleon radii are very 
similar in all cases except for ^He. The VMC wave functions underbind by about 4 (6) 
MeV in the A=Q (7) nuclei compared to GFMC; all indications are that this deficiency is 
due to a small amount of high-energy (~1 GeV) contamination, probably due to inadequate 
short-range many-body correlations. In turn, the GFMC energies are about 4% short of 
experiment, indicating the AV18/UIX model is not quite satisfactory for these nuclei. 

The contributions of the different components of weak vector and axial- vector charge and 
current operators to the reduced matrix elements (RME's) contributing to the /3-decay of 
^He and the e-capture in ''Be are reported in Tables || and |VII| . In these tables the column 
labeled "1-body" lists the contributions associated with the one-body terms of the charge 
and current operators, including relativistic corrections proportional to These are 

the operators given in Eqs. ([4. 41), (|4.7D, ([4. 9]), and ( |4.1(j| ) of Sec. |V[ The column labeled 
"Mesonic" lists the contributions from two-body vector and axial- vector charge and current 
operators, associated with pion and vector- meson exchanges, namely those of Eqs. (4.16)- 
(4.17), (4.30)-(4.31), (4.32)-(4.34), and (4.35)-(4.37) of Ref. |T3|. Lastly, the column labeled 
"A" lists the contributions arising from A excitation, obtained in perturbation theory and 
in the static A approximation, as in Eqs. (4.44), (4.48), (4.50) and (4.52) of Ref. [0. We 
reiterate that the coupling constant g\ in the A^A axial two-body current has been set 
equal to the value 1.17 qa, required to reproduce the tritium Gamow- Teller (GT) matrix 
element in a calculation based on on the same treatment of A degrees of freedom, and 
using correlated-hyperspherical-harmonics trinucleon wave functions corresponding to the 
AV18/UIX Hamiltonian model (see discussion in Sec. fV\) . 

The contributions of the different components of the weak axial current to the GT matrix 
elements occurring in the ^He /9-decay and ^Be e-capture are reported in Tables |ITT| , 0, 
and VI. The notation in these tables is similar to that just discussed above, with the only 
difference that it is now in reference to the axial current only. Furthermore, the one-body 
contributions are separated into the contributions associated with the leading and next-to- 
leading terms in the non- relativistic expansion of the covariant single-nucleon axial current, 
Eqs. (liAlp and (^121) . 

Having clarified the notation in (most of) the tables, we now proceed to discuss the 
results for the ^He /3-decay and ''Be e-capture separately in the next two subsections. 



A. The ^He /3-decay 

The ^He /3-decay is induced by the weak axial- vector charge and current, and weak 
vector current operators via the multipoles Ci{q; A), Li{q; A), Ei{q; A), and Mi(g; V). The 
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values for the associated RME's, obtained with Type I VMC wave functions, are reported 
in Table |l|at a value 0.015 fm~^ of the lepton momentum transfer q. Note that in the decay 
q varies, ignoring tiny recoil corrections, between and ~ 0.020 fm~^, corresponding to an 
end-point energy m(^He)-m(^Li)=4.013 MeV/c^. 

The largest (in magnitude) RME's are the Li{A) and Ei{A) due to transitions induced by 
the weak axial current j(q; A). This is to be expected since the decay ^He — *• ^Li e~ is a su- 
perallowed one, having {Ji'; Tj) = (0'^;l) and T/) = (1"'";0). The two-body "Mesonic" con- 
tributions to j(q; A) are individually very small, moreover they interferfere destructively and 
nearly cancel out. The two-body axial current contributions due to A excitation are at the 
level of ~ 1.5% of the leading one-body contributions. 

The transitions induced by the axial charge and vector current are inhibited, since the 
first non-vanishing terms in the long- wavelength expansions of the associated multipole 
operators are linear in g, and hence their contributions Ci{A) and MiiV) suppressed by 
one power of qR {R is the nuclear radius) with respect to Li{A) and Ei{A). The two-body 
contributions to Mi{V) and Ci{A) are relatively large, and increase (in magnitude) the one- 
body results by 9% and 25%, respectively. Among the "Mesonic" terms, the vr axial charge 
and vector current operators are dominant, while the vector-meson exchange as well as A 
operators play a minor role. It is important to stress the model-independent character of 
the TT-exchange two-body operators, whose vector and axial-vector structures are dictated, 
respectively, by gauge invariance [jl3| and current algebra [T^ arguments. 



The Mi{V) RME is approximately related by CVC to the RME of the electromagnetic 
multipole operator Mi (7) connecting the ^Li*(3.56 MeV) excited state with {.P]T) assign- 
ments (0"'";1) to the ^Li(g.s.) ground state ( J'^; T) = (1"'";0) via 

Mi(g; V) = -^('Li, 10| [t± , y ■ jt (gx; 7) 1 fRe, 00) 

~-v^Mi(g;7), (6.1) 

where we have made use of Eqs. ([4.3|) and ( ^.4|) , and have assumed that the ^Li(g.s.) as well 
as the ^He(g.s.) and ^Li*(3.56 MeV) states are members, respectively, of an isosinglet and an 
isotriplet. Of course, electromagnetic terms and isospin-symmetry-breaking components in 
the strong- interaction sector, both of which are present in the AV18/UIX Hamiltonian model, 
will spoil the relation above to some extent. However, the present VMC wave functions of 
Type I for ^Li(g.s.), ''Li*(3.56 MeV), and 6He(g.s.) are pure T = and TMt=10 and 1-1, 
respectively. Similar (Type I) wave functions have been recently employed in Ref. |11| to 



carry out a calculation of the elastic and transition form factors of the A=6 systems and, 
in particular, of the radiative width of the ^Li*(3.56 MeV) state. The values for the Mi (7) 
RME were found to be —12.81 x 10^'^ and —13.09 x 10^^ including one-body only and both 
one- and two-body terms in the electromagnetic current operator [|lll] in agreement, on the 



basis of Eq. ( |6.1| ), with those reported in Table Incidentally, the radiative width of the 
^Li*(3.56 MeV) was predicted to be lTl| 7.49 eV and 9.06 eV with one- and (one+two)-body 



currents. The experimental value is (8.19 it 0.17) eV. 

In Table |T| we list results for the GT RME, related to the Li{q = 0; A) and Ei{q = 0; A) 
RME's via 

E,{q = 0;A) = V2L^{q = 0;A) = -^g^GT , (6.2) 

V 67r 



15 



with (7^=1-2654. A few comments are in order. Firstly, the predicted ^He GT RME is 



about 5% larger than the value of 2.173 that is derived from Eqs. ( p.5|) and (|5.14|) using 
the most recent tabulation of the log(/t) value for the ^He decay, 2.910 ± 0.002, reported in 
Ref. [^. This over-prediction is already present at the level of the one-body contributions; 



those associated with two-body operators further increase the discrepancy from about 3% to 
5%. Secondly, the difference between the results obtained with Type I and II wave functions 
are very small, even though those for ^He have very different nucleon radii. The Type II 
wave functions, in contrast to those of Type I, incorporate long-range Coulomb correlation 
effects and the correct two-body clustering behavior in the asymptotic region, as discussed 
in Sec. 0. However, these asymptotically improved wave functions have only a marginal 
impact on the value of the GT RME, by reducing it by only ~ 0.2%. Thirdly, the relativistic 
corrections to j(q; A) (proportional to 1/m^) are relatively large, comparable to the leading 
two-body contributions associated with A excitation, and have been neglected in all previous 
studies we are aware of. Lastly, the Ei{q = 0; A) RME, derived from Eq. ( |6.2| ), is —i 0.6657 
(total), which should be compared to the value —i 0.6654 (total. Type I) obtained at g=0.015 
fm~^. Thus retardation corrections are tiny, as expected. 



Finally, in Table |^ we list the values for the half-life of He derived from Eqs. ( ^.5|) 
and ( [5. 61 ) under different approximation schemes for purpose of illustration (note that 
ri/2=ln2/r). The first row in Table |I^ is obtained by retaining only the Li{A) and Ei{A) 
RME's evaluated at q=0, namely neglecting retardation corrections as well as the contri- 
butions from transitions induced by the vector current and axial charge. It is equivalent 
to using Eqs. (|5.14|) and ( ^.15|) , apart from negligible recoiling corrections. The second 
row again includes only the Li{A) and Ei{A) RME's, but now keeps their full momentum 
transfer dependence. The third row includes all contributing RME's with their intrinsic 
g-dependence. The last row reports the measured half-life from Ref. [BH . 

It is important to stress that in the present calculation the effects of Coulomb distortion of 
the outgoing electron wave function are considered within the approximate scheme discussed 
in Sec. 0. In existing tabulations, such as those in Ref. these effects are treated more 
realistically, by solving the Dirac equation for the electron in the field generated by the 
(extended) charge distribution of the daughter nucleus. Obviously, this approach complicates 
considerably the formulae derived in Sec. |I| using plane waves. However, the "error" made in 
our present treatment should not be large, as can be inferred from the following argument. 
If one takes the experimental value for the GT matrix element (2.173) to compute back the 
experimental half-life, but using Eqs. ( p. 5] ) and ( p. 6] ), one obtains a value of 820.4 ms, which 
should be compared to the "true" experimental half-life of 806.7 ms. This 1.6% difference 
is presumably due to our present approximate scheme for dealing with Coulomb distortions 
of the electron waves. Lastly, if these were to be altogether ignored [by setting F{Z,pJ=l 
in Eq. ( |5.6|) ], the resulting calculated values for Ti/2, in the same approximation as in the 
third row of Table |rV|, would be 823.9 ms and 803.1 ms with 1- and (l-|-2)-body operators, 
respectively. 



B. The ^Be e-capture 

The ''Be nucleus decays by electron capture to the ground state of ^Li and to its first- 
excited state at 0.478 MeV. The (J'';T) assignments of ^Be(g.s.), ^Li(g.s.), and ^Li*(0.48 
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MeV) are (| ; i), (| ; |), and ; i), respectively. 

The calculated GT matrix elements for the transitions to the ground and first excited 
states of ^Li are given in Tables and The first and second rows in Table ( |VI|) 



list the results obtained with VMC wave functions of Type I, and based on random walks 
consisting of 20,000 (15,500) and 10,000 (10,000) configurations, respectively. Note that the 
statistical Monte Carlo errors are at the 1% level for the decay to ^Li, and almost an order of 
magnitude smaller for the decay to ^Li*. Indeed, the central values computed in the longer 
and shorter random walks are fully consistent (within errors) in the ground- to excited-state 
transition, but barely so in the ground- to ground-state transition. The third row in Tables 



and ^ lists the results obtained with VMC wave functions of Type II and based on 10,000 
point random walks. As for the A=6 case, it appears that the asymptotically improved wave 
functions of Type II lead to values for the GT matrix elements not statistically different from 
those of Type I. The two-body contributions increase by roughly 3% the one-body matrix 
elements, and as discussed below, bring theory into better agreement with experiment. 
The e-capture to ^Li(g.s.) also proceeds through a Fermi- type transition. The Fermi 



matrix element, defined as in Eq. ( p.l6| ), is F=—J2Jf + 1 with wave functions of Type I, 



in which isospin-symmetry-breaking components are ignored. However, the Type II wave 
functions do include long-range Coulomb correlations, and therefore break the isospin sym- 
metry of the isodoublet ''Be(g.s.)-^Li(g.s.). The Fermi matrix element is calculated to be, in 
this case, -1.999. 



In Table [VI I| we report the contributions from the individual components of the weak 
vector and axial charge and current operators to the dominant RME's occurring in the 
transition to the first-excited state of ''Li (the corresponding neutrino energy is E^=0.384: 
MeV or Ei,=0.752 in units of the electron mass). The Mi{V) and Ci{A) transition strengths 
are down by 3 and 4 order of magnitudes with respect to the leading Ei{A) and Li{A). 
There are in principle additional contributions from order 2 multipoles, such as C2{V) and 
M2{A), for example, however, these are expected to be even more suppressed than those 
due to Mi{V) and Ci{A). No attempt has been made to calculate them. One should note 
that the retardation corrections in the Ei[q]A) and Li{q] A) RME's are negligible. Indeed, 
Ei{q=0; A)=^/2 Li{q=0; A)=6A56 x 10"^ in the (l+2)-body calculation, which should be 
compared with 6.427 x 10~^ from Table |VI1] . Lastly, the Mi{V) RME can be expressed via 
CVC, again ignoring isospin-symmetry-breaking effects, as 

Mi(g; V) ^ -i^£.[CBe* \\ /ii(7) iT Be) - (^Li* || \\' Li)] , (6.3) 

where '^Be* is the first excited state of '^Be at 0.429 MeV with (J'^; T) = (l-/2; 1/2), and 
yLii(7) is the magnetic moment operator. From the experimentally known radiative widths 
of the ^Be* and ^Li* states ||3^, the isovector combination above of transition magnetic 



moments is found to be (5.87 ±0.14) jj,^, and the resulting Mi{q; V) is — i 5.66 x 10~^, which 
should be compared to the predicted values of — i 5.010 x 10~^ and — i 5.904 x 10~^ with one- 
and (one+two)-body currents from Table |V11| , respectively. Incidentally, in the case of the 
transition to the ground state of ''Li, one finds, using the experimental values for the ^Be 
and ^Li magnetic moments (respectively, -1.398(15) /i^v and 3.256424(2) /i^v from Ref. [§§]), 
that the isovector combination similar to that in Eq. ( |6.3|) (with ^Be* — >^Be and ^Li* — »''Li) 
is i (0.1007± 0.0003) x 10"^ while the calculated values are i 0.111 x 10^^ and i 0.132 x 10"^ 
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with one- and (one+two)-body currents. 

Finally, in Table [VIII] we report the predicted half-life and branching ratio of ''Be 
to the ground and first-excited states of "^Li. We have ignored the tiny corrections due 
to retardation effects and transitions induced by the vector current and axial charge 
operators, i.e. Eq. ( |5.19| ) has been used for F", but for the replacement |-Ris(0)p — *• 
Bx |_Ris(0)p=7.2403x 10~^. The half-life is over-predicted by about 9%, while the branching 
ratio is under-predicted by 1%. Two-body contributions reduce significantly the discrepancy 
between the calculated and measured quantities. 



VII. CONCLUSIONS 

In the present study we have reported on calculations of the ^He /?-decay and ^Be e- 
capture rates, using variational Monte Carlo (VMC) wave functions derived from a realistic 
Hamiltonian, and a realistic model for the nuclear weak current and charge operators, con- 
sisting of one- and two-body terms. Both processes are superallowed, and are therefore 
driven almost entirely by the axial current (and, additionally, by the vector charge in the 
case of the ''Be decay to the ground state of ^Li). The two-body part in the axial current 
operator has been adjusted to reproduce the experimentally known Gamow- Teller (GT) 
matrix element in /3-decay. 

The GT matrix element in ^He is over-predicted by about 5%, while those in ''Be connect- 
ing to the ground- and first-excited states of ^Li are under-predicted by 5%, when compared 
to the experimental values. However, the observed branching ratio in the ^Be e-capture is 
reasonably well reproduced by theory. We have verified explicitly that the (relatively small) 
discrepancies between the measured and calculated GT matrix elements are not explained 
by the inclusion of retardation effects or by the shift of strength to suppressed transitions 
induced by the weak vector current and axial charge operators. 

One- and two-body axial current contributions interfere constructively, leading to a 
1.7% (4.4% on average) increase in the one-body prediction for the GT matrix element 
in ^He (^Be). As a result, the inclusion of two-body operators has the effect of slightly 
worsening (significantly improving) the agreement between theory and experiment in ^He 
(^Be) systems. It is important to stress that the same model for the nuclear weak current 
adopted here, has been recently shown to provide an excellent description of the process 
^B.e{fi~ , u^^H [17|, in which two-body (vector and axial-vector) operators contribute 12% 



of the total rate. 

The origin of the current unsatisfactory situation between theory and experiment is 
most likely to be in the approximate character of the VMC wave functions used here. We 
have explored the sensitivity of the results to alternative VMC wave functions, denoted as 



Type II in Sec. |T|, which incorporate a better treatment of the asymptotic behavior, in 
particular the clustering properties into 2+4 or 3+4 sub-systems. No significant changes 
in the calculated values have been found. Thus, the next logical step in our quest for a 
quantitative understanding of weak transitions in the A=6 and 7 (as well as A=8) systems, 
is to repeat the present calculations with the more accurate Green's function Monte Carlo 
wave functions P |T(]|J5B| . 

Another, although probably lesser problem, is the inadequacy of the AV18/UIX Hamil- 
tonian, as reflected in its underbinding of the A=6 and 7 nuclei. Thus it will also be 
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worthwhile repeating these calculations with one of the more advanced Illinois three-nucleon 
potentials . It will also be useful to investigate the numerical implications of models for 
the nuclear weak current derived from effective field theory approaches, such as those in 
Refs. 



37,38 



Most of the computational techniques and computer codes developed here can 
be carried over to this planned next stage. 

Finally, it is important to establish whether the present approach, based on realistic 
interactions and currents, leads to a consistent description of the available experimental 
data on weak transitions in light nuclei, beyond the ^4=3 systems for which its validity has 
already been ascertained. One important implication of such a program should be made 
clear: it would serve to corroborate the robustness of our recent predictions for the cross 
sections of the proton weak captures on |1l2|,^ and ^He ||13|J38[ . 
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TABLES 



TABLE I. Binding energies (MeV) and point radii (fm) of A=6 and 7 nuclei obtained with 
VMC wave functions of Type I and II for the AV18/UIX Hamiltonian model. Also listed are the 
corresponding GFMC and experimental values. 



Nucleus 


Wave function 


B 


(r2)i/2 


(r2)V2 

\ 71/ 


^He(g.s.) 


Type I 




23.99(7) 


1.97(1) 


2.88(2) 




Type II 




23.78(7) 


2.28(3) 


3.23(5) 




GFMC 




28.1(1) 


1.97(1) 


2.94(1) 




Expt 




29.27 






^Lifg.s.) 


Type I 




27.09(7) 


2.49(1) 


2.49(1) 




Type II 




27.34(7) 


2.50(1) 


2.50(1) 




GFMC 




31.1(1) 


2.57(1) 


2.57(1) 




Expt 




31.99 


2.43 




'''Befg.s.) 


Type I 




30.49(9) 


2.44(1) 


2.30(1) 




Tvnp TT 




30.00(10) 


2.44(1) 


2.33(1) 




GFMC 




36.2(1) 


2.52(1) 


2.33(1) 




Expt 




37.60 






'Li(g.s.) 


Type I 




32.09(9) 


2.30(1) 


2.44(1) 




Type II 




31.87(9) 


2.32(1) 


2.43(1) 




GFMC 




37.8(1) 


2.33(1) 


2.52(1) 




Expt 




39.24 


2.27 




^Li*(0.48 MeV) 


Type I 




31.97(9) 


2.29(1) 


2.44(1) 




Type II 




31.43(10) 


2.31(1) 


2.42(1) 




GFMC 




37.5(2) 








Expt 




38.76 






TABLE II. 


Contributions to the (purely imaginary) reduced matrix elements 


EiiA), L,iA), 


Mi{V), and Ci{A) at g=0.015 fm-^ 


in ^He /?-decay. See text for notation. 




RME 


1-body 


Mesonic 


A 


(l-F2)-body 


10^ X Ei{A) 


-6.540(15) 


-0.004(3) 


-0.110(3) 


-6.654(16) 


10^ X Li(A) 


-4.623(11) 


-0.003(2) 


-0.078(2) 


-4.704(11) 


10^ X Mi(F) 


3.922(9) 


0.300(2) 


0.057(1) 


4.279(10) 


10^ X Ci{A) 


-4.584(32) 


-0.977(10) 


-0.171(5) 


-5.733(34) 


TABLE III. 


Contributions to the Gamow-Teller matrix element in ^He /?-decay. See text for 


notation. 












Wave function 


1-body NR 


1-body RC 


Mesonic 


A 


(l+2)-body 


Type I 


2.254(5) 


-0.0094(3) 


0.0014(9) 


0.0376(10) 


2.284(5) 


Type II 


2.246(10) 


-0.0100(3) 


0.0011(10) 


0.0418(11) 


2.278(10) 
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TABLE IV. Values in ms for the ^He half-life obtained in a number of approximation schemes. 
See text for an explanation. The measured half-life is also listed. 





1-body 




(l-h2)-body 


Ei{0;A) and Li{0;A) 
Ei{q;A) and Li{q;A) 
All 
Expt 


762.9 
763.5 
764.4 




743.2 
744.0 
745.1 
806.7 ±1.5 


TABLE V. Contributions to the Gamow- Teller matrix element in the ^Be e-capture to the ^Li 
ground state. See text for notation. 


Wave function 1-body NR 


1-body RC Mesonic 


A 


(l-h2)-body 


Type I(20k) 2.366(29) 
Type I(lOk) 2.288(42) 
Type Il(lOk) 2.321(41) 


-0.038(2) 0.0039(18) 
-0.034(3) -0.0029(26) 
-0.041(2) -0.0008(9) 


0.110(3) 
0.110(5) 
0.108(5) 


2.441(29) 
2.361(42) 
2.387(41) 


TABLE VI. Contributions to the Gamow- Teller matrix element in the ^Be e-capture to the 
^Li first excited state. See text for notation. 


\^'Mve funclicni l-lxKiy XR 


l-l)ody RC Mesouic 


A 


(l+2)-body 


Type I(15.5k) 2.157(6) 
Type I(lOk) 2.156(7) 
Type Il(lOk) 2.154(7) 


0.028(1) 0.0064(8) 
-0.028(1) 0.0063(10) 
-0.034(1) 0.0049(11) 


0.080(1) 
0.080(2) 
0.093(2) 


2.215(6) 
2.215(7) 
2.218(8) 


TABLE VII. Contributions to the (purely imaginary) reduced matrix elements Ei{A), Li{A), 
Mi{V), and Ci{A) in the ^Be e-capture to the ^Li first excited state at £^,^=0.7515 (in units of the 
electron mass). See text for notation. 


RME 1-body 
10^ X Ei{A) 6.192(18) 
10^ X Li{A) 4.405(13) 
10^ X Mi{V) -5.010(13) 

lO'"^ X C\(A) 5.328(90) 


Mesonic 
0.016(3) 
0.017(2) 
-0.746(6) 

0.-15(3(20) 


A 

0.220(4) 
0.163(3) 
-0.148(2) 

0.103(17) 


(l-F2)-body 
6.427(19) 
4.584(13) 
-5.904(14) 

5.918(93) 


TABLE VIII. The half-life and branching ratio of ''Be to the ground and first excited states of 
^Li, predicted with one- and (one-l-two)-body currents, are compared with the experimental values. 
The VMC wave functions of Type I are used. 


1-body 


(l+2)-body 


Expt 


Ti/2(days) 62.29 
C (%) 10.20 


58.24 
10.33 




53.22 ± 0.06 
10.44 ± 0.04 
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